Two-dimensional expansion of a condensed dense Bose gas 
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We study the expansion dynamics of a condensate in a strongly interacting Bose gas in the presence 
of an obstacle. Our focus is on the generation of shock waves after the Bose gas has passed the 
obstacle. The strongly interacting Bose gas is described in the slave-boson representation. A saddle- 
point approximation provides a nonlinear equation of motion for the macroscopic wave function, 
analogous to the Gross-Pitaevskii equation of a weakly interacting Bose gas but with different 
nonlinearity. We compare the results with the Gross-Pitaevskii dynamics of a weakly interacting 
Bose gas and find a similar behavior with a slower behavior of the strongly interacting system. 
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I. INTRODUCTION 



Many interesting features of a trapped Bose gas arc 
described by the Gross-Pitaevskii equation. This in- 
cludes dynamical studies of the Bose-Einstein condensate 
(BEG), such as the expansion of a BEG after switching 
off the trapping potential or the interference of two sep- 
arate BEG's, has been performed with remarkable accu- 
racy [l|]. The Gross-Pitaevskii (GP) equation provides 
the macroscopic quantum state of the condensed part of 
the Bose gas. It can be derived from a microscopic statis- 
tical model of many-body states at temperature T in the 
limit T = by a saddle-point approximation. However, 
the GP equation is restricted to a dilute Bose gas, with 
less then one particle in the scattering volume. More re- 
cent experiments with optical lattices have revealed that 
a much richer physics appears in a dense (or strongly in- 
teracting) Bose gas . The idea is that a static periodic 
potential of the optical lattice is provided by counter- 
propagating Laser fields, where particles occupy the lo- 
cal minima of the periodic potential. As soon as there is 
more than one particle per minimum the Bose gas must 
be considered as strongly interacting. An immediate ef- 
fect of stronger interaction is the depletion of the con- 
densate caused by the collision of particles. This is also 
known from the famous example of interacting bosons 
in form of superfluid helium. If the density of the Bose 
gas increases further we can even destroy the conden- 
sate completely and create a new quantum state in the 
form of a Mott insulator. In contrast to a BEG, the 
Mott insulator is characterized by local conservation of 
the particle number (i.e. n = 1,2,... particles per mini- 
mum of the optical lattice) but without phase coherence. 
In a trapped Bose gas, contrary to a translational in- 
variant Bose gas, both states can co-exist in the same 
system, which is known as the wedding-cake structure: 
at commensurate densities the system is in a Mott state 
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with constant density and at incommensurate densities 
the system is in a Bose-Einstein condensed state with 
spatially changing density. 

These strongly interacting systems cannot be described 
within the GP approach because the latter only takes 
into account the condensate and neglects the interaction 
with non-condensed particles. In particular, depletion 
of the condensate or the formation of a Mott insulating 
state is not accessible by the GP approach. Fortunately, 
there is an extension of the GP equation that is able to 
take into account the interaction with the non-condensed 
particles. It is known as the slave-boson (SB) approach 
and provides also a nonlinear Schrodinger equation for 
the macroscopic quantum state of the condensed part of 
the Bose gas. The interaction with the non-condensed 
particles leads to a modified nonlinearity though. For a 
dilute Bose gas the nonlinearity is the same as that of 
the GP equation but in the regime of the dense Bose gas 
it is weaker than that of the GP equation. The equation 
of the strongly interacting Bose gas also indicates the 
disappearance of the condensate when we approach the 
Mott-insulating state at higher density. Previous studies 
of a dense trapped Bose gas have shown that the BEG 
can be completely destroyed at the trapping center due 
to depletion at higher densities [3,|j]. On the other hand, 
the vortex structure is not much affected in the strongly 
interacting system because the density of the condensate 
is low in the vicinity of the vortex core [7| . 

The two-dimensional expansion of a dilute BEG past 
an obstacle is a subject of intensive theoretical and ex- 
perimental studies. A relevant parameter in these stud- 
ies is the Mach number M. It is defined as the ratio of 
the asymptotic velocity of the flow to the sound velocity 
in the medium [8|. For subsonic Mach numbers in the 
range 0.5 ^ M < 0.9, it was reported the generation of 
pairs of vortices and antivortices [^ . In case of a super- 
sonic flow, the generation of oblique dark solitons inside, 
and Kelvin ship waves outside the Mach cone (imaginary 
lines drawn from the obstacle at angles ±arcsin(l/M) 
with the horizontal axis) were found |1CHL4| . It was 
shown in [l^ that these dark solitons are convectively 
unstable for large enough flow velocity (M > 1.5), 
i.e., practically stable in the region around the obstacle. 
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FIG. 1: Behavior of the nonlinear term — x in Eq. ([2]) as a 
function of |<i>| for /j. — and ^ — 10. 



A general theory on dispersive shock waves for super- 
sonic flow past an extended obstacle was developed in 
[16| and a review paper on dark solitons in BEC can be 
found at [13 ■ Experiments addressing this problem were 
described in p8l - |20 |. Recently, a renewed theoretical in- 
terest in this issue was brought by the observation of an 
alternating vortex emission for a suitable set of param- 
eters. These are analogous to the "von Karman vortex 
street" in classical dissipative fluids j2]| . 

In this paper we shall study the effect of the interac- 
tion between the BEC and the non-condensed particles 
in a dynamical situation, where a BEC is released from 
a parabolic trap and passes an obstacle. The obstacle is 
modeled by an impenetrable disc. Due to a complex in- 
terference the macroscopic wave function will experience 
strong density fluctuations. The results of our numerical 
simulation, based on the strongly interacting gas (SIG) 
equation, will be compared with previous calculations, 
based on the GP equation. The paper is organized as 
follows: We briefly introduce the SIG equation and com- 
pare it with GP equation. Then we present the results of 
our numerical simulation for an expanding cloud in two 
dimensions that passes an obstacle. Finally, we discuss 
these results and compare them with those of the GP 
approach. 



II. SLAVE-BOSON APPROACH 

We start from a Bose gas with hard-core interaction of 
a given radius a. Then the Bose gas can be approximated 
by a lattice gas with lattice constant a. In other words, 
a provides the shortest relevant length scale in our Bose 
gas. This scale remains in the Bose gas even after its 
release from the trap. In general, a strongly interacting 
Bose gas has two constituents, namely condensed parti- 
cles and non-condensed particles. This is the case even 
at zero temperature. The interplay of all these parti- 
cles can be described by the slave-boson approach 043 ■ 
Although this is a many-body picture, the macroscopic 
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FIG. 2: Numerical solution of Eq. ([6]) for the initial profile 
given in Eq. ((7| with ai = 16 and (T2 — VSai. The density 
plot is shown for i = 2 and fi = 0. We observe the generation 
of pairs of vortices and antivortices past the obstacle. An 
impenetrable disc of radius r = 1 is placed at (25, 0) in a 
BEC radially expanding from the centre (0, 0). 



wave function of the condensate is extracted by a varia- 
tional procedure with respect to the density of the con- 
densate analogous to the derivation of the GP equation 
from the weakly interacting many-body Bose gas model. 
The corresponding effective Hamiltonian of the macro- 
scopic condensate state $(r) is 
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J (ai - a2x) , 



(1) 



where A is the three-dimensional Laplacian. Moreover, 
X is a term that is nonlinear in |<&p: 
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where Z{r) is the integral expression 

„ 2 sinh7(r) 



Z{r) = 



7(r) 



dip 



with 



(r) = ^{^ + i,l2f + |$(r) 



7(r 



(2) 



(3) 



(4) 



The coefficients are ai = f + a, a2 = cx{\ -I- \j o?\ where 
a is a numerical constant a ~ 1/5.5 (3|, |7|. The kinetic 
parameter J is associated with the mass of the bosons m 
by the relation 
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The dimensionless parameter /.t is a one-particle chemical 
potential that is associated with the density of bosons. 
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FIG. 3: Numerical solution of Eq. ((G]) for the initial profile 
given in Eq. with gi = a2 = 8. The density plot is shown 
for f = 2 and /i = 0. The dark "V" structure corresponds to a 
pair of oblique dark solitons and the oscillation in front of the 
obstacle corresponds to the "ship" waves. An impenetrable 
disc of radius r = 1 is placed at (25, 0) in a BEG radially 
expanding from the centre (0, 0). 
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FIG. 4: Cross sections of the density distributions for t = 2 
at a:: = 60 for chemical potential ^ = —1, ^. = 0, ji = 1 and 
fj, = 2. The profile for different small n is almost the same. 



This can be understood as a potential that controls the 
exchange of particles with the gas outside the trapped 
cloud by assuming that the cloud is a grand-canonical 
ensemble of atoms. Then the value of ^ fixes the number 
of bosonic atoms in equilibrium. 

The obstacle is included by choosing specific boundary 
conditions for the macroscopic wave function. In our 
case this is a hard disc, where the wavef unction vanishes 
inside the disc. Finally, the number of condensed bosons 
A^o is determined by an integral of |$(r)p over the entire 
volume of a three-dimensional Bose gas as 



a-^ (1 + 1/a)^ J 



(5) 



A. Discussion of the nonlinear term 

The form of the nonlinear term x is crucial for the 
physics of the strongly interacting Bose gas and makes all 
the difference in comparison to the GP approach. Despite 
the fact that it is defined in the integral representation of 
Eqs. (I2|), ([3]), and @, the interpretation of its asymptotic 
behavior is relatively simple. First of all, this term can 
be expanded for low condensate density (i.e. for small $) 
to get 

x~x(o) + x'(o)|i>|' 

with x'(0) < Q- The increasing behavior of this term 
with |(f>p reflects the repulsive nature of the interaction 
between the bosons. The truncation after the quadratic 
term |<I>p leaves us with the Gross-Pitaevskii Hamilto- 
nian. In the opposite limit, i.e. for |<I>| '^ oo, we get 
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The behavior of — % as a function of |$p is shown in 
Fig. [T] for two different values of the chemical potential 
yU. The nonlinearity is obviously much weaker than that 
of the GP equation. Moreover, it becomes weaker with 
increasing /i (i.e. increasing density of the Bose gas). 
This reflects the depletion of the condensate for increas- 
ing densities. 



B. Dynamics of the macroscopic wave function 

The dynamics of the quantum state of the condensate 
is given by the nonlinear Schrodinger equation: 



~6~ 



A$(r,t) 4- J(ai-a2X)*(r,t), (6) 



which will be called SIG equation. This differential equa- 
tion describes the expansion of the condensed part of 
a bosonic cloud for a given initial state $(r,0) at time 
t = 0. The spatial scale is given by scattering length 
Us ^ a. A typical experimental value for ®^Rb atoms 
near a Feshbach resonance is as ~ 200 nm Q . 



III. NUMERICAL SOLUTION AND RESULTS 

To study the evolution of an initial state $(r, i = 0) 
with Eq. (|n]) we assume a two-dimensional (2D) situa- 
tion, where the expansion is possible in the x, y direction 
but not in the z direction. This is a typical case in which 
the trapping potential, after the formation of the BEG, 
is switched off in two directions but a strongly confining 
potential is kept in the perpendicular direction. Then we 
choose for the initial state a 2D Gaussian function by 



$(a;, y,t = 0) = y^ cxp -0.5 
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FIG. 5: Numerical solution of Eq. ([6]) for the initial profile 
given in Eq. with a\ = (T2 =8. The density plot is shown 
for i = 6 and ^ = 10. The dark "V" structure corresponds to a 
pair of oblique dark solitons and the oscillation in front of the 
obstacle corresponds to the "ship" waves. An impenetrable 
disc of radius r = 1 is placed at (25, 0) in a BEG radially 
expanding from the centre (0, 0). 
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FIG. 6: Numerical solution of Eq. ([6]) for the initial profile 
given in Eq. ((7|) with cri = (72 = 8. The density plot is shown 
for f = 6 and ^ = 20. The dark "V" structure corresponds to a 
pair of obhque dark solitons and the oscillation in front of the 
obstacle corresponds to the "ship" waves. An impenetrable 
disc of radius r = 1 is placed at (25, 0) in a BEG radially 
expanding from the centre (0,0). We see that for bigger ^l, 
the expansion of the BEG is much slower than in the case 
illustrated in Fig. [S] 



where tri and (T2 are parameters which determine the 
shape of the Gaussian. 

We apply a 2D finite-difference method (Crank- 
Nicolson method) combined with a split-step method to 
solve equation ^ numerically. Schematically this reads 



$(x,2/,t) = e 

g-iV^-^-iTSt 



-iV^ 



-^^*$(x,2/,0) 
f$(a;,y,0) 

5tl 



(8) 
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where V and T are the potential and kinetic term, re- 
spectively, 5t is the time step and CJ\f\x, 5t] is the action 
of the Crank-Nicolson method for the kinetic term on the 
X direction. This approach is accurate up to 0{5t^), and 
good results arc obtained for sufficiently small time step 



5t. We use a Gauss-Hcrmite quadrature to calculate at 
each time step the integrals x and Z{r) in ^ and (O, 
respectively, and noticed that the Gauss-Hermite inte- 
gration converges already for only nine points. 

In our simulations, we use as boundary condition 
<I>(r) = at the boundary of a square region. An im- 
penetrable obstacle is implemented by assuming that the 
order parameter vanishes inside a circle of radius r = 1. 
The obstacle can also be implemented as a strong repul- 
sive Gaussian potential ^(r) = C/oexp(— (.x^ + y'^)/2n) 
with the chemical potential given by /^(r) = ^i — V{r). 
Both conditions give similar results. 

It was shown in [0] that the flow of a dilute BEG past 
an obstacle generates vortices and antivortices at sub- 
sonic velocity. Increasing the flow velocity, the frequency 
of vortices generation also increases. For large enough 
flow velocity, the vortices generation is so often, that the 
distance between them becomes less than their radial size 
and it takes a long time for their separation from each 
other at finite distance from the obstacle |10l - KL3 , Il5| . We 
study this transition from generation of vortices to gen- 
eration of dark solitons in a dense Bose gas, considering 
different profiles of the initial state (different values of cri 
and (72 in Eq. (O), which corresponds to faster expan- 
sion for a small initial profile and to slower expansion, 
for bigger initial profile. 

We observe the emergence of pairs of vortices and an- 
tivortices using (Ti = 16 and cr2 = \/3cri in Eq. ^ (cf. 
Fig. [5]). This is similar to what was observed in recent 
experiments [20|. 

We study the case of faster expansion of the dense 
BEG for the initial state in Eq. ([7]) with cti = (T2 = 8 
and different values of /i. namely ^ = — l,/i = 0,/i=l, 
/i = 2, /^ = 10 and ji = 20. We observe in Fig. ([3] - 
ini) a pair of oblique dark solitons pas t the obstacle and 
linear waves, as in the case of GP ^M,- For small values 
of chemical potential, /i = —1, /i = 0, /.t = 1, /i = 2, the 
pattern (position, slope and amplitude of the soliton) is 
almost the same (cf. Fig. S]). However, we find that the 
velocity of the expansion of the cloud becomes smaller 
for an increasing chemical potential (see Fig. [5] and Fig. 

We also study the effect of the nonlinear term x in 
Eq. (O with the full potential of Eq. ^ for different 
values of /i on the dynamics. Since x becomes smaller as 
we increase ^ (see Fig. [71), the expansion velocity of the 
cloud is also reduced with increasing ^. Gorrespondly, 
the effective potential term of the Hamiltonian is strongly 
repulsive for small ^ and becomes flat for bigger /i (see 
Fig. El). 

Although the nonlinearity of the SIG equation is more 
complex than its counterpart in the GP equation, the 
flow past an obstacle shows qualitatively the same pat- 
terns as for the GP equation: it is characterized by 
oblique dark solitons and by linear waves. Our study 
suggests that these effects shall well appear in a strongly 
interacting gas too and may be observed experimentally. 




IV. CONCLUSION 



FIG. 7: Cross section at a:: = of the nonlinear term x in Eq. 
((2]) for ^ = and fi — 10. It goes to zero as we increase fj.. 
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FIG. 8: Cross section at a; = of the potential term in Eq. 
(|6]) for n — and /i — 10. It becomes flat as we increase fi. It 
makes the expansion of the cloud slower for big /i. 



We have studied the two-dimensional expansion of a 
strongly interacting BEC in the presence of an obsta- 
cle. The strong interaction was treated within the slave- 
boson approach, which leads to a nonlinear Schrodingcr 
equation with a special nonlinearity, different from the 
Gross-Pitacvskii term. We solved numerically this equa- 
tion and analyzed the contribution of the nonlinearity. 
Similar to what was found in previous studies of the GP 
equation [3, Il3l , we observed the pairwise generation of 
vortices and antivortices (at low velocity) and of shock 
waves (at high velocity), respectively. We noticed that 
bigger values of the chemical potential (i.e. for higher 
densities of the Bose gas) the expansion is slower due to 
a reduced nonlinear term. The characteristic features of 
the SIG equation in terms of a generation of vortices, 
oblique dark solitons and ship waves are similar to those 
described by the GP. 
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